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A largely unsolved theoretical issue in controlled fusion research is the consistent kinetic treatment 
of slowly-time varying plasma states occurring in coUisionless and magnetized axisymmetric plasmas. 
The phenomenology may include finite pressure anisotropics as well as strong toroidal and poloidal 
differential rotation, characteristic of Tokamak plasmas. Despite the fact that physical phenomena 
occurring in fusion plasmas depend fundamentally on the microscopic particle phase-space dynamics, 
their consistent kinetic treatment remains still essentially unchalleged to date. The goal of this paper 
is to address the problem within the framework of Vlasov-Maxwell description. The gyrokinetic 
treatment of charged particles dynamics is adopted for the construction of asymptotic solutions 
for the quasi-stationary species kinetic distribution functions. These are expressed in terms of 
the particle exact and adiabatic invariants. The theory relies on a perturbative approach, which 
permits to construct asymptotic analytical solutions of the Vlasov-Maxwell system. In this way, 
both diamagnetic and energy corrections are included consistently into the theory. In particular, by 
imposing suitable kinetic constraints, the existence of generalized bi-Maxwellian asymptotic kinetic 
equilibria is pointed out. The theory applies for toroidal rotation velocity of the order of the ion 
thermal speed. These solutions satisfy identically also the constraints imposed by the Maxwell 
equations, i.e. quasi-neutrality and Ampere's law. As a result, it is shown that, in the presence 
of non-uniform fluid and EM flelds, these kinetic equilibria can sustain simultaneously toroidal 
differential rotation, quasi-stationary finite poloidal flows and temperature anisotropy. 



PACS numbers: 52.55.Fa, 52.30.Gz, 52.25.Xz, 52.25.Dg 

I. INTRODUCTION 

Plasma dynamics is most frequently treated in the 
framework of stand-alone MHD approaches, i.e., formu- 
lated independent of an underlying kinetic theory. How- 
ever, these treatments can provide at most a partial de- 
scription of plasma phenomenology. The reason is re- 
lated to two basic inconsistencies of customary fluid ap- 
proaches. First, the set of fluid equations may not be 
closed, requiring in principle the prescription of arbitrary 
higher-order fluid fields. Second, in these approaches typ- 
ically no account is given of microscopic phase-space par- 
ticle dynamics as well as phase-space plasma collective 
phenomena. It is well known that only in the context of 
kinetic theory these difficulties can be consistently met. 
Such a treatment in fact permits to obtain well-defined 
constitutive equations for the relevant fluid fields describ- 
ing the plasma state, overcoming at the same time the 
closure problem. Kinetic theory is appropriate, for exam- 
ple, in the case of coUisionless or weakly-coUisional plas- 
mas where phase-space particle dynamics is expected to 
play a dominant role. 

Unfortunately, for a wide range of physical effects aris- 
ing in magnetically-confined plasmas and relevant for 
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controlled fusion research, a fully consistent approach 
of this type is still missing. Surprisingly, these include 
even the description of equilibrium or slowly-time vary- 
ing phenomena occurring in realistic laboratory Tokamak 
plasmas. The issue concerns specifically the description 
of finite pressure anisotropics, strong toroidal differential 
rotation as well as concurrent poloidal flows observed in 
Tokamak devices. The deficiency may represent a serious 
obstacle for meaningful developments in plasma physics 
(both theoretical and computational) and controlled fu- 
sion research. In particular, it is well-known that both 
toroidal and poloidal plasma equilibrium rotation flows 
may exist in Tokamak plasmas ^L, 2j. The observation 
of intrinsic rotation, occurring without any external mo- 
mentum source @, remains essentially unexplained to 
date, being mostly ascribed to turbulence or boundary- 
layer phenomena occurring in the outer regions of the 
plasma [I,!!]. Such an effect, potentially combining both 
toroidal and poloidal flow velocities with temperature 
anisotropy, may be of critical importance both for sta- 
bility and suppression of turbulence ^-S^. 

The goal of the present investigation is the con- 
struction of slowly-time varying particular solutions of 
the Vlasov-Maxwell system for coUisionless axisymmetric 
plasmas immersed in strong magnetic and electric fields. 
In principle, two approaches are possible for the investi- 
gation of the problem. One is based on the Chapman- 
Enskog solution of the drift-kinetic Vlasov equation, 
namely achieved by seeking a perturbative solution of the 
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form fs = Jms + e/is + where < e ^ 1 is an appro- 
priate dimensionless parameter to be defined below (see 
Section 2) and /ms a suitable equilibrium kinetic distri- 
bution function (KDF). In customary formulations this is 
typically identified with a drifted Maxwellian KDF. An 
example is provided by Hinton et at. Q where an ap- 
proximate equilibrium KDF carrying both toroidal and 
poloidal fiows was introduced to describe ion poloidal 
flows in Tokamaks near the plasma edge. An alterna- 
tive approach is represented by the construction of ex- 
act or asymptotic solutions of the Vlasov equations of 
the form fs = f*s, with /*s to be considered only func- 
tion of particle exact and adiabatic invariants, via the 
introduction of suitable kinetic constraints. This tech- 
nique is exemplified by Ref.liol], where /*s was assumed 
to be a function of only two invariants, namely the par- 
ticle energy = ^^6$*^ and toroidal canonical momen- 
tum Pips = ^^4'*s [see their definitions given below], and 
identified with a generalized Maxwellian distribution of 
the form 

Here H^g is the invariant H^s = Es — J d'ipilo{''P) i 



while A.^,s = {n*s_T*s} denotes suitable "structure func- 
tions", i.e., properly defined functions of the particle in- 
variants. In Refs.[l0l-[l2l| these were prescribed imposing 
the kinetic constraint A*s = A*s('(/'*s)- By performing a 
perturbative expansion in the canonical momentum (see 
also the related discussion in Section 6), it was shown 
that /*s recovers the Chapman-Enskog form, with the 
leading-order Maxwellian KDF carrying isotropic tem- 
perature Ts (?/;), species-independent toroidal angular ro- 
tation velocity ^o{ip) (see definition given by Eg. ((39)) ) 
and finite toroidal differential rotation, i.e., -^Qo{ip) ^ 
0. A basic aspect of Tokamak plasmas is the property of 
allowing toroidal rotation velocities Rflo comparable to 

1/2 

the ion thermal velocity ut/ii = {2Ti/Mi} ' . As shown in 
Ref.^io'l this implies the fundamental consequence that, 
for kinetic equilibria characterized by purely toroidal dif- 
ferential rotation as described by the KDF ([T]), neces- 
sarily the self-generated electrostatic (ES) potential $ in 
the plasma must satisfy the ordering (Mivff^^) / (ZiC^) ^ 
0{e). If the ion and electron temperatures are compa- 
rable, in the sense that Ti/T^ ~ O (e°), it follows that 
an analogous ordering must hold also for the electron 
species. Therefore, the same asymptotic condition must 
he adopted for all thermal particles of the plasma, namely 
for which |v| ~ Vtus, independent of species. 

In the following, utilizing such a type of ordering, the 
second route is adopted. Hence, the theory developed 
here applies to a two-species ion-electron plasma charac- 
terized by toroidal rotation velocity of the order of the ion 
thermal speed. It relies on the perturbative kinetic the- 
ory developed in Refs.[l^,[l^ (hereafter referred to as Pa- 
pers I and II). The aim is to provide a systematic gener- 



alization of the theory presented in Ref.H^, allowing /^^ 
to depend on the complete set of independent adiabatic 
invariants, and therefore to vary slowly in time ("equilib- 
rium" KDF). In particular, here we intend to show that, 
besides the properties indicated above, also temperature 
anisotropy, finite poloidal flow velocities and flrst-order 
perturbative corrections, including finite Larmor-radius 
(FLR) corrections, can be consistently dealt with at the 
equilibrium level. A remarkable feature of the approach 
is that, by construction, all the moment equations stem- 
ming from the Vlasov equation are identically satisfied, 
together with their related solubility conditions (i.e., fol- 
lowing from the condition of periodicity of the KDF and 
its moments in the poloidal angle). An interesting devel- 
opment consists in the inclusion of both diamagnetic (i.e., 
FLR) and energy corrections arising from the Taylor- 
expansions of the relevant structure functions. In such 
a case the structure functions are identified with smooth 
functions of both the particle energy and toroidal canon- 
ical momentum, of the general form 

Ks^ As{i'*s,^*s). (2) 

with the functions As(-0, $) being identified with suitable 
fiuid fields, s denoting the species index. This permits 
the construction of a systematic perturbative expansion 
also for the KDF itself, allowing to retain perturbative 
corrections (of arbitrary order) expressed as polynomial 
functions in terms of the particle velocity. In particu- 
lar, under suitable assumptions, the leading-order KDF 
is shown to be determined by a bi-Maxwellian distribu- 
tion carrying anisotropic temperature and non-uniform, 
both toroidal and poloidal, fiow velocities. Thanks to the 
kinetic constraints, constitutive equations are determined 
for the related equilibrium fiuid fields. First-order correc- 
tions with respect to e are shown to be linear functions 
of suitably-generalized thermodynamic forces. These in- 
clude now, besides the customary ones [l^, additional 
thermodynamic forces associated to energy derivatives 
of the relevant structure functions. 

The constraints imposed by the Maxwell equations are 
then investigated. First, the Poisson equation is an- 
alyzed within the quasi-neutrality approximation. As 
a development with respect to Ref.|lC|, it is proved 
that the perturbative scheme determines uniquely, cor- 
rect through O (e°) , the equilibrium ES potential, in- 
cluding the I/O (e) contribution. Second, the solubility 
conditions of Ampere's law are shown to prescribe con- 
straints on the species poloidal and toroidal flow veloc- 
ities and the corresponding current densities. The the- 
ory applies for magnetic configurations with nested and 
closed toroidal magnetic surfaces characterized by finite 
aspect ratio. 

The paper is organized as follows. First, in Sec- 
tions 2 and 3 the Vlasov-Maxwell and magnetized-plasma 
asymptotic orderings are posed, together with the basic 
assumptions concerning the plasma and its electromag- 
netic (EM) field. In Section 4 particle first integrals and 
adiabatic invariants are recalled, including guiding-center 
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adiabatic invariants predicted by gyrokinetic theory. In 
Section 5 particular solutions of the coUisionless Vlasov 
equation are investigated. Their representation in terms 
of suitable structure functions is discussed. Then, based 
on the Taylor-expansion of the relevant structure func- 
tions, in Section 6 a perturbative kinetic theory is ob- 
tained for the KDF. As an application, the leading-order 
and the first-order diamagnetic and energy contributions 
to the KDF arc displayed. In Section 7 the connection 
between the kinetic and fluid treatments is addressed. In 
Section 8 the leading-order number density and flow ve- 
locity carried by the stationary KDF are reported. The 
implications of Maxwell equations are discussed in Sec- 
tions 9 and 10. In particular, in Section 9 the issue con- 
cerning the quasi-neutrality condition is addressed. It is 
shown that quasi-neutrality determines uniquely, up to 
an arbitrary constant, the ES potential (THM.l) and is 
consistent with the plasma asymptotic orderings intro- 
duced (Corollary to THM.l). Then, constraints placed 
by the Ampere equation are investigated in Section 10 
(THM.2). Relevant comparisons with previous literature, 
based either on kinetic or fluid approaches, are presented 
in Section 11. Finally, concluding remarks are given in 
Section 12. 



II. VLASOV-MAXWELL ASYMPTOTIC 
ORDERINGS 

In the following, for particles belonging to the s- 
species, we introduce the characteristic time and length 
scales Ats = and ALs = AL = 27rr, with 27rr 

1 /2 

and v±ths = {T±s/Ms} denoting respectively the 
connection length and the thermal velocity associated 
to the species perpendicular temperature T±s (defined 
with respect to the local magnetic field direction). We 
shall consider phenomena occurring in time intervals Ats 
which belong to the ranges Tps <C Atg ^ tcs, where 

^ \4nn^fz~ej^) ' ^^'^ isotropic specics temper- 

atures tcs = ^^^"^ \nA(z e)^ denote respectively the 
Langmuir time and the Spitzer ion self-collision time. 
A similar ordering follows for the corresponding scale- 
length ALs letting ALg = AtgVths^ with Vths being the 
species isotropic-temperature thermal velocity. For defi- 
niteness, we shall consider here a plasma consisting of n 
species of charged particles, with n > 2. Such a plasma 
can be regarded, respectively, as: 

(#1) CoUisionless: in validity of the inequality be- 
tween Atg and tcs , contributions proportional to the ra- 
tio ecs = 1; here referred to as the collision-time 
parameter, can be ignored. Thus, Coulomb binary inter- 
actions arc negligible, so that all particle species in the 
plasma can be regarded as coUisionless. 

(#2) Continuous: thanks to the left-side inequality 
between Atg and Tps , plasma particles interact with each 
other only via a continuum mean EM field. In particular. 



the inequality SLg,s = ^ 1 is assumed to hold, with 
^Lg,s denoting the Langmuir-time parameter. 

(#3) Quasi-neutral: due again to the same inequal- 
ity, the plasma is quasi-neutral on the spatial scale ALg 
corresponding to Atg. 

Systems fulfilling requirements #l-#2 - the so-called 
Vlasov- Maxwell plasmas - rely on kinetic theory, since 
fluid MHD approaches are inapplicable in such a case 
(see related discussion in Papers I and II). Such plasmas 
are described in the framework of the Vlasov-Maxwell 
kinetic theory. In this case the plasma is treated as 
an ensemble of particle s— species (subsets of like par- 
ticles) each one described by a KDF fg{z,t) defined 
in the phase-space F = x r„ (with F^ C and 
F„ = denoting respectively the configuration and ve- 
locity spaces) and satisfying the Vlasov kinetic equation. 
Velocity moments of fs(z,t) are then defined as inte- 
grals of the form jp d^vQ{z,t)fg{z,t), with Q{z,t) be- 
ing a suitable phase-space weight function. In particular, 
for Q{z,t) = {l,v} the velocity moments determine the 
source of the EM self-field {E'*^'-^, B'*'^'-^}, identified with 
the plasma charge and current densities {p{r,t),J{r,t)}. 

In addition, we require the plasma to be axisymmet- 
ric, so that, when referred to a set of cylindrical coor- 
dinates {R,(fi,z), all relevant dynamical variables char- 
acterizing the plasma (e.g., the fluid fields and the EM 
field) are independent of the azimuthal angle ip. Here, 
by assumption, the configuration space is identified with 
the bounded internal domain of an axisymmetric torus, 
which can be parametrized in terms of the scale-lengths 
(a, Rq), with a denoting a = sup {r, r G F^} and Rq the 
radius of the plasma magnetic axis. 

III. BASIC ASSUMPTIONS 

In this section the basic hypothesis of the model, which 
include the EM field and the magnetized-plasma order- 
ings, are pointed out. 

A. The EM field 

Here we restrict our analysis to EM fields which are 

slowly-time varying in the sense [E(x, e'^t), B(x, e'^t)] , 
with A; > 1 being a suitable integer {quasi-stationarity 
condition). This type of time dependence is thought to 
arise either due to external sources or boundary condi- 
tions. In particular, the magnetic field B is assumed to 
be of the form 

B = V X A = B'^'^f (x, sH) + B«^* (x, e'^t) , (3) 

where B*^'-'' and B®^* denote the self-generated mag- 
netic field produced by the plasma and a finite external 
magnetic field produced by external coils. In particular 
the magnetic field B admits by assumption a family of 

nested and closed axisymmetric toroidal magnetic sur- 
faces {'ip{ x)} = {tp{x) = const.}, where ^ denotes the 
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poloidal magnetic flux of B and, because of axisymme- 
try, X can be identified witfi tfie coordinates x = {R, z). 
In sucli a setting a set of magnetic coordinates ("0, tp, "&) 
can be defined, wliere is a curvilinear angle-like co- 
ordinate on the magnetic surfaces V'(x) — const. It is 
assumed that the vectors (VV', Vip, Vi9) define a right- 
handed system. Each relevant physical quantity G'(x, t) 
can then be conveniently expressed either in terms of the 
cylindrical coordinates or as a function of the magnetic 
coordinates, i.e. G(x, <) = G{ip,'0,t). The total mag- 
netic field is then decomposed as 



B = /(x, eH)\7ip + Vi/'(x, e'^t) x Wip, 



(4) 



where Bt = I{x,e'''t)\/ip and Bp = V'i/;(x, e'^t) x Vip are 
the toroidal and poloidal components of the field. In par- 
ticular, the following ordering is assumed to hold: ^ 
O (e") . Finally, the corresponding electric field expressed 
in terms of the EM potentials {<I>(x, e^^i), A(x, e'^t)} is 
considered primarily electrostatic, namely 



1 dA 



E(x,£'=i) = -V$ -e'^'-— ^ -V$ 

C OT 



(5) 



with T denoting the slow-time variable r = e'^'i, and 
quasi-orthogonal to the magnetic field, in the sense that 
^ ^ 0{e), while ^ O (eO) . Together with 

the quasi-stationarity condition, this implies that, to 
leading order in e, <I> = ^{4', £^t). In particular, assuming 
that both <I> and A are analytic with respect to e, it can 
be shown that, consistent with gyrokinetic (GK) theory 
and the asymptotic orderings indicated below (see next 
Section), they must be considered of the general form 

$ = (V^,e''^t) -He°$o (V',^,e''^) + (6) 



where $ is expressed in terms of the magnetic coordinates 
and A_i is A_i = '0V(/9 -I- g (^tp, i?, e'^tj Vt9, with g being 
a suitable function. 



B. The magnetized-plasma orderings 

Next, let us introduce the magnetized plasma ordering 
appropriate for the treatment of single-particle dynam- 
ics in magnetized plasmas, i.e. for which in particular 
^ E^. For s = i,e, this requires the definition of the 
following additional dimensionless parameters: 

1) Larmor-radius parameter eM,s = 317" ^'^'^ Larmor- 
time parameter eLr,s = -^f- ■ here t^s and r^g are re- 
spectively the Larmor time and the Larmor radius of the 
species s, with s — 1, n, defined as rj^^ = v±ths/^csi with 
flcs = ZgcB / MsC = 1/tls being the species Larmor fre- 
quency. Imposing the requirement that t^s Aig and 



rLs ^ ALs, it follows that Em.s and eLr,s are infinitesi- 
mals of the same order, i.e., < eLr,s ~ £m,s ^ 1- Re- 
quiring again that Ti ^ Te, and furthermore Zi ^ O (1), 

1 /2 

it follows that eM,i ~ (if^) ^M,e- 
2) Canonical-momentum parameter: Ss 



14, 



where = v • and L^s denotes the species 



where 



particle angular momentum. 

3) Total-energy parameter: i 

-^w^ ~ Ts and ZsC^ axe respectively the particle kinetic 
and ES energy. 

In principle, the parameters and as are independent 
(in particular, as pointed out in Paper I, they might differ 
from eM,s)- More precisely, here we shall consider the 
subset of phase-space for which the following ordering 
holds: 



< a. 



£M.s ^ 1, 



(8) 



which applies in the subset of thermal particles. No- 
tice that the assumption on Eg is consistent with the re- 
quirement of finite inverse aspect-ratio (see below), while, 
as recalled above, the ordering on ag is required for the 
treatment of Tokamak equilibria in the presence of strong 
toroidal differential rotation [lol - [l3j . The same orderings 
are of course invoked also for the validity of the GK the- 
ory (see Ref.[l3] and also Eg. dTSl) in the next section 
and the related discussion). The assumption on the ag- 
ordering can be shown to be consistent with the quasi- 
neutrality condition (see Corollary to THM.l in Section 
9). The previous requirements imply, for all species, the 
asymptotic perturbative expansions in the variables ■0*s 
and $*,: 



= -0 [1 + O (Em,*)] : 

= $[l + 0(eM,s)] 



(9) 
(10) 



Finally, to warrant the validity of the Vlasov equa- 
tion on the Larmor-radius scale, we shall impose also 
that < £mfp,s ~ £Cs < £m,s, with Smfp.s = ^ and 
£cs = denoting respectively the mean-free-path pa- 
rameter and the collision-time parameter. Then, consis- 
tent with quasi-neutrality, we demand also that £Lg,s ~ 
£d < £m,s < e < 1, with e — sup {ea/.s, s = e,i} . Fi- 
nally, the inverse aspect-ratio parameter 6 = will be 

considered finite, i.e. such that 5 ^ O (e^) . We remark 
that the parameters {cts, e^, e^^.s, ea/.s} deal with the 
single-particle dynamics, {eLg,s,£D,s,£mfp,s,£Cs} con- 
cern collective properties of the plasma, while 5 is a 
purely geometrical quantity. 



IV. THE PARTICLE ADIABATIC INVARIANTS 

For single-particle dynamics, the exact first integrals 
of motion and the relevant adiabatic invariants are well- 
known. In particular, the adiabatic invariants can be 
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defined either in the context of Hamiltonian dynamics 
or GK theory [T7l - [l9j . In both cases, for a magnetized 
plasma, they can be referred to the Larmor frequency. 
Hence, by definition, a phase-function Ps depending on 
the s-species particle state is denoted as adiabatic in- 
variant of order n with respect to sai,s if it is con- 
served asymptotically, namely in the sense jp- In Pg — 

+ 0(eM s where n > is a suitable integer and fl'^^ 
is the Larmor frequency evaluated at the guiding-center 
position x'. Note that, in the following, we shall use a 
prime " ' " to denote a dynamical variable defined at the 
guiding- center position r' (or x' in axisymmetry) . Under 
the assumptions of axisymmetry, the only first integral 
of motion is the canonical momentum p^ps conjugate to 
the azimuthal angle ip: 



Z,e 



Z.e 



p^s = MgR-v ■ -\ V = V'*s- (11) 



-'u:> 

c c 
Furthermore, the total particle energy 



M, 



-V 



-Zge$(x,£"t) = Z,,e$, 



(12) 



with n > 1, is assumed to be an adiabatic invariant of 
order n. 

Let us now analyze the adiabatic invariants predicted 
by GK theory. As usual, the GK treatment involves the 
construction - in terms of an asymptotic perturbative ex- 
pansion determined by means of a power series in Em.s - 
of a difFcomorphism of the form z = (r, v)— > z' = (r', v'), 
referred to as the GK transformation. The GK trans- 
formation is performed on all phase-space variables z = 
(r,v), except for the azimuthal angle tp which is left un- 
changed and is therefore to be considered as one of the 
GK variables. Here, by definition, the transformed vari- 
ables z' [GK state) are constructed so that their time 
derivatives to the relevant order in ea/.s have at least one 
ignorable coordinate, to be identified with a suitably- 
defined gyrophase 0'. Starting point is then the repre- 
sentation of the particle Lagrangian in terms of the hy- 
brid variables z. This is expressed as >Cs(z, -^z,e'^t) = 
r • P, - -H,.(z,e'=i), where P, = [Af,v + ^A{x,eH)] 
and 'Hs{'z,e'^t) = ^v'^+Zse^{x,e'^t) denotes the corre- 
sponding Hamiltonian function in hybrid variables. The 
development of GK theory is well known. It involves a 
phase-space transformation to a local reference frame in 
which the particle guiding-center is instantaneously at 
rest with respect to the ip-swciace to which it belongs. 
In this case, the leading-order GK transformation can be 
proved to be necessarily of the form 



This coincides with the so-called frozen-in velocity, 
namely the fluid velocity with respect to which each 
line of force is carried into itself. The rest of the nota- 
tion is standard. Thus, u' and w' denote respectively 
the parallel and perpendicular (guiding-center) veloci- 
ties, with w' = w' cos(p'e'i + w' smip'e'2 and 4>' denoting 
the gyrophase angle, fl'^^ = ^fj^ and b' b(x',£'^i), 
with b(x, e''t)= B(x, £'^t)/i3(x, e'^'i). Notice that, here, 
by construction, 
with respect to 



must be considered of O {sm,s) 

while for thermal particles and 
|w'| are all of the same order of Vths- In particular, due 
to the previous orderings, for the validity of GK theory 
the EM potentials (<&, A) entering the Lagrangian must 
be considered of the form indicated above (see Eqs.® 
and ([7])), namely both of 1/0 (e) with respect to the re- 
maining terms. As a consequence, the ordering (jS]) for Ug 
necessarily applies, under the assumption Ti/T^ ^ O (e'^) 
considered here. On the other hand, as in Ref . pLOl| . |U'| 
is to be taken of the order of the ion thermal velocity 
Vfhi, while |U'| ~ S^qR, with being the toroidal an- 
gular rotation frequency, defined below by Eq. ((39|) . It is 
important to stress here that these two conditions imply 
that $ must satisfy the asymptotic ordering given above 
by Eq.(E^. Therefore, the previous orderings for as and 
(f> must be regarded as basic prerequisites for the descrip- 
tion of Tokamak plasmas characterized by toroidal rota- 
tion speeds comparable to the ion thermal velocity. 

By construction, in GK description the gyrophase an- 
gle is ignorable, so that the magnetic moment is an 
adiabatic invariant of prescribed accuracy. In particular, 
the leading-order approximation is — "' ~ ^^^^ 



Ms 



2B' 



Two further adiabatic invariants can immediately be ob- 
tained from the previous considerations. In fact, since 
the azimuthal angle ip is ignorable also in GK theory, 
the conjugate GK canonical momentum p'^^, referred to 
as the guiding- center canonical momentum, is necessarily 
an adiabatic invariant. Neglecting corrections of 0(eM,s) 
this is given by 



B' 



u'l' 



+ —i>', (15) 



which provides a third-order adiabatic invariant. We re- 
mark that both m'g and p'^^ can in principle be identified 
with adiabatic invariants of 0(£^,/^ ) , with k >1 arbitrar- 
ily prescribed [2^. In the following we shall make use of 
the local invariants (V'+sj Es,m'g) to represent the particle 
state, while adopting p'^^ to deal with the dependences 
in terms of u' . 



V = u'h' 



U', 



(13) 



Here, in particular, U' = U(x',e''i), with U(x, e'^i) being 
the fluid-field identified with the E x B— drift velocity: 



U(x,e'=t) 



B 



V$ X b. 



(14) 



V. VLASOV KINETIC THEORY: 
EQUILIBRIUM KDF 

Let us now proceed constructing asymptotic solutions 
of the Vlasov equation holding for coUisionless Tokamak 
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plasmas in validity of the previous assumptions. The 
treatment is based on Papers I and II, where equilib- 
rium generalized bi-Maxwellian solutions for the KDF 
were proved to hold for accretion disk plasmas. In par- 
ticular, the following features are required for the equi- 
hbrium KDF: 

1) For all of the species, different parallel and 
perpendicular temperatures are allowed (temperature 
anisotropy) . 

2) Non-vanishing species dependent differential 
toroidal and poloidal rotation velocities are included. 

3) The KDF is required to be an adiabatic invariant 
asymptotically "close" to a local bi-Maxwellian. Hence, 
in particular, in the case of a locally non-rotating plasma 
(i.e., for which both toroidal and poloidal rotation ve- 
locities vanish identically on a given ^/'-surface) the KDF 
must be close to a locally non-rotating bi-Maxwellian. 

In analogy with Papers I and II, it is possible to show 
that Requirements 1) - 3) can be fulfilled by a suitable 
modified bi-Maxwellian expressed solely in terms of first 
integrals and adiabatic invariants, including also a suit- 
able set of structure functions {A*^} of the form ^ (see 
the precise definition below). Hence, the desired KDF is 
identified with an adiabatic invariant of the form 



(16) 



with n > 1, and where the brackets (■(/'*s, <l'*s) de- 
note the dependence in terms of the structure functions 
{A*s($*s, ■!/'*s)}- In particular, in agreement with as- 
sumptions 1) - 3), f^s is identified with KDF of the form: 



/3* 



exp 



Til 



(17) 

which we refer here to as the generalized bi-Maxwellian 
KDF with parallel velocity perturbations. The notation 
is as follows. First, {A*s} = {/3*s, a*s, Ty^g, ^^*s, are 
structure functions subject to kinetic constraints of the 
type ([2]), assumed analytic functions of both ijj*s and $*s- 
These are, by definition, suitably close to appropriate 
fluid fields A., = As{ip, $). In particular, the functions A^ 



are defined as {A^} = = ^^,( 



B' 
At. 



Til 



where r/g denotes the pseudo-density, T^^g and T±s the 
parallel and perpendicular temperatures, with = 

jT-, while rtg and are the toroidal and parallel 

rotation frequencies. Second, the phase-function E^,s is 
defined as i?*., = ^^*s ^ p'^sC*si while H^.^ is identified 
with 



H-^^s — E^ ^>^sr2^ 

c 



(18) 



We stress that the form of /*s [see Eg. dTTl) ] is ob- 
tained consistent with assumption 3), namely such that 
when the constraint fl^g = = locally holds, fi,s re- 
duces to the non-rotating generalized bi-Maxwellian KDF 



y/2 exp|-^ -m'.a^sj. In par- 



ticular, unlike Ref.[10|. the definition given above for H^, 
follows by requiring that i?*s, and hence also iJ*s, is a 
local linear function of the frequencies J7*s and and 
of the canonical momenta p^s and p' . 

An equivalent representation for (|17l) can be obtained 
invoking the previous definitions. This yields: 



/3*s exp 



T|u, 



{2T,/Msf'^ (Til,,) 



,1/2 



(19) 



X exp ■ 



(- 



2T|L 



where W*s ~ e^Rfl^s, U!,_^ = -i?^*^ and 



.7/;17„-Z,e$ + T:,, (20) 



TL = 1 



2 y 

AfscVV' • V'$' Zse 



+ ) ^*s. (21) 



Note that U'^^g is non-zero only if the toroidal magnetic 
field is non- vanishing. 

The following comments are in order: 

1) /*s is by construction a solution of the asymptotic 
Vlasov equation 



1 d 
dt 



In/,, =0-fO (£"+!). 



(22) 



2) fts is defined in the phase-space F = Fr x Fu, where 
Tr and Tu are both identified with suitable subsets of the 
Euclidean space M^. In particular, is non-zero in the 
subset of phase-space where the adiabatic invariants p'^g , 
H'g and are defined. It follows that is suitable for 
describing both circulating and trapped particles. 

3) The velocity moments of /*s, to be identified with 
the corresponding fluid fields, are unique once f^s is pre- 
scribed in terms of the structure functions. 



VI. PERTURBATION THEORY 

In this section we develop a perturbative kinetic the- 
ory for the KDF f^^. This is obtained by performing 
on a double- Taylor expansion for the implicit func- 
tional dependences in the variables "0*8 and $*s carried 
only by the structure functions {A*,}, while leaving un- 
changed all the remaining phase-space dependences. As 
indicated above, such asymptotic expansions can be ex- 
pressed in terms of the dimensionless parameters as and 
Eg in validity of the ordering ([5]) . Hence the double Taylor 
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expansion yields: 

A„ ^ As + iip*s - ^) 



and 



dA, 



9$. 



(23) 



where both and the partial derivatives in ([23]) are by 
construction functions depending only on (ip,^). This 
implies also their general dependence in terms of the mag- 
netic coordinates (ip,^) (see Section 9). We notice that 
the asymptotic order of the "gradients" of the structure 
functions and depends whether in A*^, ■i/'+s 

and/or are considered "fast" or "slow" variables with 
respect to s, in the sense that the same gradients can 
be considered respectively 0(e") or 0{e). In principle, 
different possible orderings are allowed for the pertur- 
bative expansion of Here we shall assume in par- 
ticular that the structure functions /3*s, a*s, Tiug have 



fast dependences, while f2*s,^*s have only slow ones. As 
a consequence, the set of derivatives | , f§^| and 

{ dijj"' ' } both taken here as 0(e). It follows 
that to first order in e the KDF f^s can be approximated 
as: 



f^s — fs 



^Ds 



''Ds 



(24) 



where the leading-order KDF fs does not depend on the 
gradients of Ag. Hence, all the informations about the 
gradients of the structure functions appear only through 
the first-order (in e) perturbations h^j^l and /i^^. These 
are denoted respectively as the diamagnetic- correction 
(see Ref.[l3l) and the energy- correction (see Paper II), 
which result from the leading-order Taylor expansions 
with respect to V'+s and $*s. In particular, the following 
results apply. First, fs is expressed as 



{27T/Msf' {T\\sf^T^ 



X exp < 



( 



M, f V - W, - U[^h 

2riic 



m's^ K25) 



and hence is identified with a bi-Maxwellian KDF with 
parallel velocity perturbations (see Paper II). In Eq.((25|) 
Ws = r2si?^V(/3 and U'^^ ~ -gy^s are related to the 
leading-order toroidal and parallel flow velocities and de- 
pend on angular frequencies of the general form fls = 
r2s(V',$) and = ^s i'^J,^)- In addition, the function 
Us is defined in terms of the pseudo-density rjs as 



Us $) = r]s (ip, $) exp 



Xs 

Til. 



(26) 



Xs 



MsU, 
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— — + — ^ 



ij^s - Zse^ + T'\, (27) 



(28) 



Second, the diamagnetic and energy-correction contribu- 
tions h^j^l and h^^l are given by 



,(2) 



1 Zse 
Ms 



Zsc tpns 

^4 — — L-3s 



+ -^^nsAA{^,■e^){29) 



2Zse [ T||5 ( 
Here y^, i = 1, 5, is defined as 

H 



C5s } v%30) 



Ais + A2 



-A^s — A2S 



1 

% " 2 



fi'sAi 



Cis + C2s [ jr- ^ 2 ^ ^ Ai'sC4s 



(31) 
(32) 
(33) 



where Hg = Es — ^^V*s^^s and the following definitions 
have been introduced: Ais 



^2s — Tm, 1 ^3s 



Oil 1 ^^^s — g^, 
dlnfls A — 3a a A — d\n(s „„J — d In 3 — 

dij ' "^48 = -g^, = ana Oi^ = -55-, 02s = 

9111^113 „ _ ainOa r< — r — ^ 

g$ J '-'3s — : '--'4s — gtj, i '-'5s — gtj, • 

The outcome of the perturbative theory is as follows: 
1) The asymptotic expansion in terms of 7/1*5 and lead- 
ing to the diamagnetic-correction /i^] is formally analo- 
gous to that presented in Ref. [l^ . The Taylor expansion 
in terms of (energy expansion) is instead a novel 
feature of the present approach and leads to the energy- 



(2) 
Da 



correction h 

2) The kinetic equilibrium /*s is compatible with 
species-dependent rotational frequencies fis and ^g. No 
restriction follows from the KDF on their relative mag- 
nitudes, so that the general ordering ^ ~ O is per- 
mitted. 

3) A fundamental feature is related to the functional 
dependences imposed by the kinetic constraints on the 
structure functions. As a basic consequence, the latter 
depend both on the poloidal flux -p and the ES potential 
$. As proved below, the ES potential $ is generally a 
function of the form $ = $(x, e'^i), with x = {R,z), i.e. 
it is not simply a i/)-flux function. Hence, when expressed 
in magnetic coordinates, the structure functions become 
generally of the form A^ = As d, e'^i) . This type of 
functional dependence is expected to apply for arbitrary 
nested magnetic surfaces having finite inverse aspect ra- 
tio. On the other hand, in the case of large aspect ratio 
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{1/5 » 1), the poloidal dependences in As are expected 

to become negligible. Nevertheless, /i]^^ remains finite 
even in this case. The reason is that also in this limit the 
double Taylor expansion (f23| still applies. 

4) The coefficients Ais and Cts, i = 1,5, can be identi- 
fied with effective thermodynamic forces, containing the 
spatial variations of Ag respectively across the ^ = const. 
and $ = const, surfaces. 



VII. THE VLASOV FLUID APPROACH 

An elementary consequence concerns the fluid ap- 
proach defined in terms of the Vlasov description, i.e., 
based on the moment equations following from the 
asymptotic Vlasov kinetic equation (see Eg.ip^'). In fact, 
assuming that the KDF is identified with the adiabatic 
invariant given by Eq. (|16p , these equations are necessar- 
ily all identically satisfied in an asymptotic sense, namely 
neglecting corrections of O (e""*"^) . Furthermore, because 
/*s is by construction periodic, also the corresponding 
solubility conditions, related to the requirement of peri- 
odicity in terms of the i9-coordinate, are necessarily ful- 
filled. To prove these statements we notice that if Q{z) is 
an arbitrary weight function, identified for example with 
Q — (l, V, u^) , then the generic moment of Eq. ((22|) is: 

^ rf'«Q^/*s = + O(e"+i), (34) 

where r„ denotes the appropriate velocity space of inte- 
gration. Using the chain rule, and taking into account 
explicitly also the dependence in terms of p'^^ , this can 
be written as 

/ d\Q JLf + lu £^ h + O (e"+^) . 

{ ' dt dm' ' dt dp' ) 

(35) 

On the other hand, Eq. ([M|) can also be represented as 
I d^^[jtiQf*s]~f*s^^Q]=0 + O{e'^+'), (36) 

which recovers the usual form of the velocity-moment 
equations in terms of suitable (and uniquely defined) 
fluid fields. For Q = (l,v) one obtains, in particular, 
that the species continuity and linear momentum fluid 
equations are satisfied identically up to infinitesimals of 
O (e"+^) . Similarly, the law of conservation of the species 
total canonical momentum can be recovered by setting 
Q = ijj^s, namely 

^ d\-^^[Asf*s]=0 + O{e'^+'). (37) 

In the stationary case this implies the customary species 
angular momentum conservation law for the species an- 
gular momentum L*°* = Msi?^n*°*V*°* • Vip (see Pa- 
per II). Here the notation is standard. In particular 



the velocity moments of the KDF | n*°* , V*°* , W ,Ll°*j 

can be introduced, to be referred to as species num- 
ber density, flow velocity, tensor pressure and canoni- 
cal toroidal momentum. They are defined by the inte- 
grals Jp d^vQfi,s, where Q is now identified respectively 

with Q = |l,-|^,Ms(v- V*°*)(v- V*°*),%^V„}. 

It is worth remarking here that the velocity moments 
are unique once the KDF f^.^ [see Eq. (jl7p ] is prescribed 
in terms of the structure functions {A*s} . On the other 
hand, as a result of Eqs. (f22|) and (p4l) . it follows that 
the stationary fluid moments calculated in terms of the 
KDF /*s are identically solutions of the corresponding 
stationary fluid moment equations. 

We conclude this section pointing out that no restric- 
tions can possibly be required on the KDF and the EM 
potentials as a consequence of the validity of these mo- 
ment equations. Therefore, the only possible constraints 
on the KDF are necessarily only those arising from the 
solubility conditions of the Maxwell equations. 



VIII. CONSTITUTIVE EQUATIONS FOR 
SPECIES NUMBER DENSITY AND FLOW 
VELOCITY 

In this section we present the leading-order expres- 
sions of the species number density and flow velocities 
predicted by the kinetic equilibrium. The calculation of 
these fluid moments is required for the subsequent anal- 
ysis of the Maxwell equations. An explicit calculation 
of the moment integrals can be carried out by adopting 
the perturbative asymptotic expansion of described in 
Section 6. This also requires to perform an inverse GK 
transformation, by expressing all of the guiding-center 
quantities appearing in the equilibrium KDF in terms of 
the actual particle position, according to Eq. p^ . 

Consider flrst the evaluation of the species flow veloc- 
ity V*°*. Adopting the GK representation for the particle 
velocity, the leading-order contribution to the flow veloc- 
ity is found to be 

Ys = U+4(f^.+6)b + 

+ + - U] • (1 - bb) ,(38) 

\\s 

where U is the frozen-in velocity deflned by Eg. (1141) . 
Then, ignoring correction of 0(e), U can be approxi- 
mated as U = i?^f7oViy9- (l — bb) . Here fio is the species- 
independent and ip-flux function [lol | 

^o{^,s'^t)^c^ (39) 

and ($) ^ § B7^$ denotes the i/)-surface average, 
with >(~^ = § b'^w • Then, in terms of the relative 
toroidal frequency Arj^ = D,s — the leading-order 
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flow velocity becomes 



where 



0„ 



-is 



[An, 



(40) 



This implies that can be decomposed in terms of the 
total toroidal and poloidal rotation velocities 



Vps 



Ps 



|Vi?| 



V, • ep, 



(41) 
(42) 



where ep = j^f]'; the corresponding rotation fre- 
quencies r^Ts and flps are respectively: 

!^Ts = Qo + [Arts + is] + 

J- II « 



rip, = [Al],+e.];g^ (l-^ 



(43) 
(44) 



with i = VV' X V(/3 • V?9. 
We remark that: 

1) To leading-order in e, the poloidal flow velocity 
(1421) is non-zero only in the presence of temperature 
anisotropy. More precisely, provided 7^ 1, a non- 
vanishing Vps may arise only if Afig -1-^, ^0. Therefore, 
even if fig coincides with the frozen-in frequency Hq, ilps 
is different from zero if ^0. 

2) The effect of the contributions Ail, and is anal- 
ogous, although their physical origins are different. In 
particular Ail, represents the departure from the frozen- 
in rotation velocity flo, while is determines the parallel 
velocity perturbation in the KDF. 

3) If the frozen-in condition is invoked, namely fig = 
flo, Eq.dID]) becomes 



V, ^ noR'w^ + is 



T±s 9 I 



1 



(45) 

which takes into account both finite poloidal rotation 
and temperature anisotropy. In case of isotropic tem- 
peratures, i.e., = 1, the previous equation provides 

a purely toroidal flow given by Vg = + is) R^^}P- 
When is = this reduces to the customary result [lOl |. 
namely = QoR'^^'P- 

Finally, we report the calculation of the number den- 
sity n*°*. Neglecting again first-order diamagnetic and 
energy-correction contributions, the leading-order species 
number density is found to be 



Vs i'ip, *) exp 



Xs - Zse<t> 



(46) 



Xs 



Ms P , ,2 Ms o 

-^^^ins + isf-^u^ 



c 



Ms cVip ■ V$ Zse 



B 



B 



^s 

is 



1 I^' 



(47) 



contains the combined contribution of the kinetic ener- 
gies carried by the rotation frequencies fls, is and the 
frozen-in velocity U. 



IX. QUASI-NEUTRALITY 

In this section we investigate the implications of the 
quasi-neutrality condition following from the Poisson 
equation. Here by quasi-neutrality we mean that the 
equation 



E 



Zscn, 







is satisfied asymptotically in the sense that 
IV -El Ojel) 



(48) 



(49) 



with £pi = <C 1 denoting the Debye- length dimen- 
sionless parameter, with Api ~ Xds = TpsVths, AL ~ 
ALs — AtgVths, and n*°' the total species-number den- 
sity. We intend to show that the first two terms in 
the Laurent expansion ([6]) of $ can be determined from 
Ea. (H5)) by prescribing n*°* to leading-order in e, namely 
in terms of Ea. (|46p . In particular the following result 
holds. 

THM.l - Explicit form of the ES potential 

Let us assume that the species KDF is defined by 
Eq. and the finite aspect-ratio ordering applies. 

Then, imposing the quasi-neutrality condition in the 
case of a two-species ion-electron plasma, the following 
propositions hold: 

Til) Correct through O {e^), the ES potential satisfies 
the asymptotic implicit equation 



where 5 (V', $) is the source term given by 



(50) 



S'(V',i?,$) 



In 



ZtT]i 



Xe 

Til. 



X^ 

Til, 



(51) 
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with rjs being the species pseudo- density and the quantity 
X, = X, (V;, $) defined by Eq.\Jl^. 

TI2) // the temperatures are non-isotropic, then the 
species pseudo- density is generally of the form rjg — 
r]s {ip,'d,^). Instead, in the case of isotropic tempera- 
tures, rjs = r]s{i^,^). 

TI3) A particular solution consistent with the kinetic 
constraints is obtained letting Zirji = rj^. 

TI4) In particular, in validity of TI3, correct through 
O (e") the ES potential $ is uniquely determined by 
Eg. \50\l and is necessarily of the form where $_i 
obeys the equation 



(52) 



while $0 is obtained subtracting <i>_i from Eg. i f5Q|) . 

PROOF - Til - The proof of the first statement 
can be obtained from Eq.(|48|) by substituting for the 
species number density the leading-order solution given 
by Ea.(|46)). TI2 - The proof follows by noting that, 
in validity of the kinetic constraint on the species 
pseudo-density is such that = t]~^'^^'^)- 

the other hand, from the kinetic constraint imposed 
on a.^,s and the prescriptions that B — B (?/;, i?) and 
T\\s — T\\s [i^T^], it must be that T_ls is necessarily of 
the type T^s = T^s{i','&T^)- Therefore, the general 
dependence of the pseudo-density is also necessarily of 
the form r/s = 77^ (t/','!?, $). On the other hand, in the 
limit of isotropic temperatures Tj_s — Ty^ — Tg (-0, $) 
and a*s = 0. The functional dependence of r]s be- 
comes therefore of the type rjs = r]s {fp, TI3 - Due 
to the arbitrariness of the structure function /3*s, it fol- 
lows that /3e and /3i can always be defined in such a way 
that^^^ = 1 even when T±i ^ T±,,. In particular, 
this constraint is consistent with the requirement that 
the ES potential vanishes identically in the absence of 
toroidal and poloidal rotations. TI4 - By definition the 
Poisson equation, subject to suitable boundary condi- 
tions, must determine completely (i.e., uniquely) the ES 
potential <I>. Therefore, Eg. ([501 yields necessarily the 
complete solution, correct through O (e*^) . In particu- 
lar, by inspecting the order of magnitude of the different 
contributions in the source term Xg ('(/', i?,$), the Lau- 
rent expansion © can be introduced. In particular, by 
retaining in Xs(-0,'!?, $) only contributions of l/0{e), 
$-1 (tp) is found to obey Eq.^. Q.E.D. 

A fundamental implication of THM.l, and in particu- 
lar of the validity of Eq.®, is to assure the consistency 
of the perturbative erg-expansion as well as the orderings 
introduced in Sections 3 and 4. In fact, let us inspect the 
order of magnitude (with respect to the parameter e) of 
the r.h.s. of Eg. ((52)) . For definiteness, let us assume that 
{ft, + - {fie + Ce) - flo, requiring r„,/r,|e ~ O (e") 
and Zi^ O (e") . Due to Eq.dSll) it follows that the order 
of magnitude of $_i is <I>_i ~ On the basis of this 

conclusion, the following statement holds. 



Corollary to THM.l - Consistency with the CTs- 
expansion. 

Given validity of THM.l and the quasi-neutrality con- 
dition, invoking the previous assumptions it follows that 
CFi CTe^ O (e). 



PROOF - First, by assumption ai 



Zie<S> 



and 



As a consequence, due to the previ- 
ous hypotheses ct; ^ a^. Furthermore, thanks to quasi- 
neutrality, it follows that 



MsVthihvthi 



Z^e^^ 



MgVthiR 



(53) 



The order of magnitude of the two factors on the r.h.s 
follows from the asymptotic ordering for the canonical- 
momentum parameter and the requirement indicated 
above that VLqR ~ vtm (see also Ref.|10|). It is concluded 

that, since by construction O (ci) ^ 0(e), 



0{e), while -^^^ ^ 0(£*'), which manifestly implies 

the thesis. Q.E.D. 

The following further remarks are useful in order to 
gain insight in the previous results. 

1) Eg. ([50]) represents the general solution holding in 
the case of a two-species plasma characterized by tem- 
perature anisotropy, poloidal and toroidal fiow velocities. 

2) In Eq.dill) all quantities Ai^^ = {r^^, Tp} can 
be considered (to leading-order in e) as only i/j-functions, 
namely of the form aI^^ = aI^^ {tp). Therefore, Eg. ((52)) 
provides an ODE for $_i {tp). 



X. THE AMPERE EQUATION 

Let us now investigate the constraints imposed by the 
Ampere law on the leading-order current densities and 
equilibrium flows. Let us consider the case of a two- 
species plasma. The following results apply. 

THM.2 - Constraints on poloidal and toroidal 
flows. 

Given validity of the asymptotic Vlasov kinetic equa- 
tion V2'^] for the species KDF defined by Eq. {11^ , 
the quasi-neutrality condition and the magnetized- 

plasma asymptotic orderings (see Section 3), for a two- 
species plasma the following propositions hold: 

T2i) The poloidal flow velocity Vps{ip,^,'^) may be 
either species- dependent or independent. In the first case 
necessarily the constraint condition 



d_ 



Y^ZsCUsVps (V',^?,$) 



= 



(54) 



must be fulfilled. In the second case, if Eq. is not sat- 
isfied, the corresponding total equilibrium current density 
must vanish identically. 
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T22) In both cases, the toroidal flow velocity remains 
species- dependent, so that the corresponding current den- 
sity is generally non-vanishing. 

T23) Both poloidal and toroidal magnetic fields can be 
self-generated by the plasma. 

PROOF - T2i - Let us consider first the component 
of Ampere's equation along the directions orthogonal to 
\/(p. This yields the following set of two scalar equations 
for the toroidal magnetic field IVip: 

1^ = —y2z,en,Vp,{i;,^,<f)[l + 0{e)], (55) 
ayj c ^ — ' 

^ = + Oie), (56) 

implying manifestly the solubility condition (1541) . There- 
fore, either the total poloidal current density is a ijj- 
function, or the poloidal flow velocity Vps $) must 

be species-independent. The first condition can always be 
satisfied by suitably selecting the species pseudo-density. 
In fact, even in validity of TI3, the species pseudo- 
density can be defined in such a way to satisfy the con- 
straint ((M)) . Therefore, excluding the null solution, a 
non-vanishing current density must appear when Vps is 
species-dependent . 

T22 - The proof of the second statement follows by 
noting that, in validity of proposition T2i, the quantity 
^jp [A^ls + £,s] may still remain species-dependent. As a 

consequence, by direct inspection of Eq. (|4T|) . it follows 
that the toroidal current density is generally non-null. 

T22 - Thanks to the previous propositions, it follows 
that both the toroidal and poloidal magnetic fields can 
be self generated. In particular, the self toroidal field re- 
quires necessarily the presence of temperature anisotropy, 
while the poloidal self field may arise even in the case of 
isotropic temperature, due to deviations from the frozen- 
in condition fls = and/or parallel velocity perturba- 
tions associated to ^s- Q.E.D. 

We briefly mention the case of a multi-species plasma. 
In fact, in collisionless systems plasma sub-species can 
be introduced, simply based on the topology of their 
phase-space trajectories. For example, different species 
can be identified distinguishing between circulating and 
magnetically-trapped particles. These components can 
in principle be characterized by KDFs carrying dif- 
ferent structure functions, and in particular different 
poloidal flow velocities. In this case both the poloidal 
and toroidal flow velocities remain generally species- 
dependent. Therefore, the corresponding current den- 
sities may be expected to be non- vanishing. 



XI. COMPARISONS WITH LITERATURE 

An interesting issue is related to comparisons with the 
literature. For what concerns the kinetic formulation, 
the relevant benchmark is represented by Ref. 10], where 



the theory of collisional transport in toroidally rotating 
plasmas was investigated. Although the conceptual foun- 
dations of the perturbative kinetic approach here adopted 
have already been exhaustively detailed in Sections 2 
to 10, it is worth analyzing some differences arising be- 
tween the two approaches. In detail, besides the inclusion 
of temperature anisotropy, parallel and toroidal velocity 
perturbations as well as the prescription of the kinetic 
constraints, the main differences with Ref. are as fol- 
lows: 

1) The first one lies in the choice of equilibrium KDF. 
This is related, in particular, to the different defini- 
tion adopted here for the dynamical variable i?*^ (see 
Ea. ([T5)) ). The motivation for this definition have been 
detailed in Section 5. Such a choice permits to ob- 
tain an explicit analytical solution for the leading-order 
ES potential $_i {ip) (see THM.l), based uniquely on 
the quasi-neutrality condition rather than imposing fluid 
constraints (see Ref.[10|). 

2) In the present approach no constraints arising from 
the moment (i.e., fluid) equations are placed on the struc- 
ture functions {A*s} [see Eq.([2|)] and consequently on 
the velocity moments of the KDF f^s- In particular, in 
our case, unlike the case of collisional plasmas treated 
in Ref. 10], the general form of the equilibrium species- 
fluid velocity is merely a consequence of the form 
prescribed for the equilibrium KDF. Therefore, it can- 
not follow from imposing the validity of fluid equations, 
but only from the solubility conditions of the Maxwell 
equations. 

3) The analysis of the Ampere equation has been car- 
ried out to investigate its consequences on the toroidal 
and poloidal species-flow velocities in the presence of 
temperature anisotropy (see TIIM.2). The discussion ex- 
tends the treatment given in Ref. jioj , where only differen- 
tial toroidal flows were retained in the kinetic treatment. 

Let us now consider, for the sake of reference, also the 
case of statistical fluid approaches. Such treatments (in- 
cluding those adopting multi-fluid formulations) typically 
do not rely on kinetic closures conditions and/or include 
FLR as well as perturbative kinetic effects, such as dia- 
magnetic and energy-correction contributions. Further 
issues include: 

1) The treatment of kinetic constraints. As shown 
here, kinetic constraints are critical for the construction 
of the KDF. They allow the structure functions to retain, 
in principle, both tp (leading-order) and "d (first-order) de- 
pendences. The correct functional form of the fluid flelds, 
arising as a consequence of the kinetic constraints, may 
not be correctly retained in customary fluid treatments 
(see for example Refs.[2ll[2^). 

2) The proper inclusion of slowly time-dependent tem- 
perature and pressure anisotropics. As pointed out here, 
the functional form of the parallel and perpendicular tem- 
perature is related to microscopic conservation laws, in 
particular particle magnetic moment conservation. On 
the other hand, fluid approaches normally ignore such 
constraints. Even when kinetic closure conditions are in- 
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voked for the pressure tensor (see for example Ref . |22l| ) . 
their vahdity may become questionable if they are not 
based on consistent equilibrium solutions for the KDF. 

3) Another example-case is provided by the kinetic pre- 
scription for the expression of the number density, here 
shown to exhibit a complex dependence in terms of the 
ES potential, centrifugal potential as well as toroidal and 
parallel frequencies (for comparison see Rcf422]). 

4) Finally, the functional form of the poloidal flow ve- 
locity may differ from what can be obtained adopting a 
two- fluid approach 2l[. In particular, in our treatment 
the toroidal and parallel rotation frequencies are con- 
sidered independent of each other, so that kinetic con- 
straints need to be imposed separately on Q.s and ^s- 
Furthermore, according to the kinetic treatment, a non- 
vanishing equilibrium poloidal flow velocity can only ap- 
pear in the presence of temperature anisotropy. 



XII. CONCLUDING REMARKS 

In this paper, a theoretical formulation of quasi- 
stationary configurations for collisionless and axisymmet- 
ric Tokamak plasmas has been presented. This is based 
on a kinetic approach developed within the framework of 
the Vlasov-Maxwell description. It has been shown that 
a new type of asymptotic kinetic equilibria exists, which 
can be described in terms of generalized bi-Maxwellian 
distributions. By construction, these are expressed in 
terms of the relevant particle first integrals and adiabatic 
invariants. Such solutions permit the consistent treat- 
ment of a number of physical properties characteristic of 
collisionless plasmas. They include, in particular, differ- 
ential toroidal rotation and finite temperature anisotropy 



and poloidal flows in non-uniform multi-species Tokamak 
plasmas subject to intense quasi-stationary magnetic and 
electric fields. The existence of these solutions has been 
shown to be warranted by imposing appropriate kinetic 
constraints for the structure functions which appear in 
the species distribution functions. By construction, the 
theory assures the validity of the fluid moment equations 
associated to the Vlasov equation. In particular, the nov- 
elty of the approach lies in the explicit construction of 
asymptotic solutions for the fluid equations in terms of 
constitutive equations for the fluid fields. The approach 
is based on a perturbative asymptotic expansion of the 
equilibrium distribution function, which allows also the 
determination of diamagnetic and energy-correction con- 
tributions. The latter are found to be linearly propor- 
tional to suitable effective thermodynamic forces. Fi- 
nally, the constraints placed by the Maxwell equations 
have been investigated. As a result, the electrostatic 
potential has been determined by imposing the quasi- 
neutrality condition. Furthermore, it has been shown 
that non-trivial solutions for the toroidal and poloidal 
species rotation frequencies are allowed consistent with 
the solubility conditions arising from the Ampere law. 
The discussion presented here provides a useful back- 
ground for future investigations. 
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